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ABSTRACT 



Gas inflowing along stellar bars is often stalled at the location of circumnuclear rings, that form an effective reservoir for massive 
star formation and thus shape the central regions of galaxies. However, how exactly star formation is proceeding within these circum- 
nuclear starburst rings is subject of debate. Two main scenarios for this process have been put forward: In the first the onset of star 
formation is regulated by the total amount of gas present in the ring with star forming starting once a mass threshold has reached in a 
'random' position within the ring like 'popcorn'. In the second star fomation preferentially takes place near the locations where the 
gas enters the ring. This scenario has been dubbed 'pearls-on-a-string'. Here we combine new optical IFU data covering the full stellar 
bar with existing multi-wavelength data to study in detail the 580 pc radius circumnuclear starburst ring in the nearby spiral galaxy 
NGC 6951. Using HST archival data together with SAURON and OASIS IFU data, we derive the ages and stellar masses of star clusters 
as well as the total stellar content of the central region. Adding information on the molecular gas distribution, stellar and gaseous 
dynamics and extinction, we find that the circumnuclear ring in NGC 6951 is ~ 1-1.5 Gyr old and has been forming stars for most of 
that time. We see evidence for preferred sites of star formation within the ring, consistent with the 'pearls-on-a-string' scenario, when 
focusing on the youngest stellar populations. Due to the ring's longevity this signature is washed out when older stellar populations 
are included in the analysis. 

Key words. Galaxies: individual: NGC 6951 - Galaxies: ISM - Galaxies: stellar content 



<^ 1. Introduction 

Secular evolution comprises all agents of change in a galaxy 
\^ that are independent of the external environment of the galaxy. 
P\| Secular evolution, unlike hierarchical clustering and mergers, 
' takes place over long timescales, near-equal to the lifetime of 
j<~>) galaxies, and is important for a galaxy's evolution in the cur- 
rent epoch. A large part of secular evolution takes place due to 
asymmetries in the gravitational potential of a galaxy. All asym- 
\ metries assist in the inward and outward motion of gas across 
. ^ ' parts of the disk. One of the asymmetries that acts the fastest 
and most severely is a large scale bar: an elongated concentra- 
5_i tion of stars and/or gas centered on the nucleus, usually of sev- 
' eral kpc size. Depending on the strength of the bar, the sound 
speed in the gas, and the central mass concentration, a circumnu- 
clear ring will form near the inner Lindblad resonance at a radius 
defined by the bar (e.g. Athanassoula 1992a; Regan & TeubenI 
l2003l:lMazzuca et alJI2008tlKim et alJl2012ll . 

A ring directly hinders further inflow towards the nucleus, 
i.e. nu clear feeding (e.g. Macieiewski 2004; Garci'a-Burillo et al. 
120051) . A circumnuclear ring is in effect the gas reservoir for the 
gas moving inward. Determination of the life times of circumnu- 
clear rings is therefore of paramount importance to understand 
the influence they have on the evolution of their host galaxies. 
The high gas densities reached in circumnuclear rings naturally 
lead to star formation. The ages of the stellar population(s) in 
circumnuclear rings may be a good indicator of the overall age 



of the ring, and thus the stellar bar itself. A caveat naturally fol- 
lows; how does star formation proceed in circumnuclear rings? 
Is there a time scale that can be related to the onset of star forma- 
tion in the ring (and that should be added to determine the ring's 
age)? 

Much theoretical and modeling work has been done on 
the formation o f circumnuclear ririgs (e.g. Ath an assoula 
1992b; Englmaier ^Shlosmanl 120001: iRegan & TbubenI 120031: 
Macieiewski.2004) . Circumnuclear gas rings are closely related 
to the more generic nuclear spirals; asymmetries in the gravita- 
tional potential induce a wave perturbation in the gas. All gas has 
angular momentum and it needs to lose this momentum before it 
can move towards the center This loss happens in the dust lanes 
and nuclear spirals. They are the locations where gas is shocked 
and loses angular momentum, allowing gas to migrate to smaller 
orbits. 

As has b e en s hown by lAthanassou 1 ([T992ah : 
iRegan & TeubenI (l2003h and others, the existence of a cir- 
cumnuclear ring inside a large scale bar depends, among other 
factors, on the orbit families that make up the circumnuclear 
region. x\ orbits are elongated along the asymmetric potential 
and are the main orbits that make up the large scale bar. X2 
orbits on the other hand, are oriented perpendicular to the bar. 
Gas orbits slowly turn from one to the other in a bar potential. 
That X2 orbits are again symmetric to the large scale stellar 
bar potential, is important, because there should be no net 
gravitational torques over the orbits (i.e. no further shocks in 
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the gas). X2 orbits are therefore able to serve as repositories for 
gas. If there is a significant fraction of .51:2 orbits present in the 
circumnuclear region, gas can settle on them and build up a 
circumnuclear gas ring. 

The connection between dust lanes, nuclear spirals and cir- 
cumnuclear rings is very clear observationally. The dust/gas 
lanes of the bar do not extend to the nucleus in the case of a ring 
and actually end at the ring radius. Two points can often be seen 
in the gaseous ring, offset by 180° , where the lanes con nect to 
the rin g. They are called 'contact points', following Bok er et al.l 
(120081). So, after being shocked and having lost most angular 
momentum, infalling gas enters the circumnuclear ring at the 
contact points and will stay on the X2 orbits. The immediate re- 
sult of this situation is that the gas density in the ring will in- 
crease with continued inflow, with star formation as the natural 
consequence. 

The most straight-forward star formation model for circum- 
nuclear rings has been proposed by Elmegreen (1994). The gas 
density will build up until it reaches a critical density for star for- 
mation. The entire ring will become unstable and star formation 
will proceed throughout the ring. The time scale of the onset of 
star formation is then directly related to the inflow rate onto the 
ring. This scenario was called 'popcorn' by B oker et al.. (2008.) . 
who themselves proposed another scenario. 

This second scenario proposes an azimuthal age gradient for 
the recently formed stars, starting at the contact points. Since gas 
enters the ring at two distinct positions, it stands to reason that 
the gas density in the ring is higher at the contact points (termed 
'overdensity regions', or 'ODR' by the authors) than elsewhere 
in the ring. The critical density for star formation will be reached 
first at (or close to) those positions. By overall rotation of the 
system these newly formed star clusters will be azimuthally dis- 
tributed in the ring, while star formation continues at the contact 
points, like 'pearls-on-a-string'. 

NGC 695 1 is a good candidate to investigate which of these 
two scenarios holds true and estimate the life tim e of the cir- 
cumnuclear ring. NGC 695 1 is a nearby (24. 1 Mpc, iTuUvll 19881) 
SAB(rs)bc galaxy. It has a circumnuclear ringat a radius of 5' 



or 580pc. It has a low incHnation (46.2° . iHaan et all l2009l) . 
which provides a clear view of the entire ring. Its molecu- 
lar gas dynamics in and around the ring have already been 
studie d in detail by some of the au thors ("Van der Laan et al.' 



201 Ih. as well as by ofliers (iKohno et al. 1999; Krips et al. 2007; 



Storchi-Bergmann et al. 2007). The gravitational torques acting 
on the cold molecular gas were also analyzed (, H aan et al. 2009: 
IVan der Laan et al.1 l201lh . iMazzuca et"al] (l2008h have investi- 
gated the Hot regions in the ring, and found young stellar ages 
(<10Myr) in the northern segment of the ring and older ages 
(~ IQMyr for an instantaneous, or ~ 100 Myr for a continuous 
star formation scenario) in the rest of the circumnuclearring. No 
other determination of the ages of the stellar populations of the 
ring has been published to date. 

New observations and archival data retrieved for this study 
will be presented in §2. In §3 the known relevant properties of 
the ring are summarized. In §4 the extinction in the ring is de- 
rived, before detailing the detection of star clusters in §5. §6 
contains the determination of the ages of the star clusters and 
determination of the ages of the overall stellar content based on 
SAURON data. A discussion with respect to star formation sce- 
narios and ring life time is given in §7. A summary follows in 
§8. 




Fig. 1. Overlay of the Field-of-View of various datasets pre- 
sented in this paper From small to large: OASIS (white), HST 
(green), PdBI (blue), SAURON (red). The underlying (B band) im- 
age was extracted from DSS. 



2. Observations and Archival Data 



2.1. OASIS data 



The OASIS Integral-Field Unit (iMcDermid et al.l l2004 is in- 
stalled at the William Herschel Telescope (WHT) in La Palma 
(Canary Islands). Two individual OASIS configurations were 
used, namely the MR 5 16 and MR 661, which, respectively, 
cover the spectral range around the H/3, Mg lines ([4760 - 
5560A]) and around Ha, [Nil], [SII] lines ([6210 -701 OA]) at a 
resolution of R ~ 1500, and a sampling of about 2. 15A per pixel 
(see the OASIS/WHT page for further details). For both config- 
urations, the 22 mm enlarger was used, providing a field of view 
of approximately KX'S x 7'.'4 (see Fig.lTJ with original spaxels of 
^.'26. 

NGC 695 1 was targeted during a run in July 2006, and we 
obtained five and four 30min exposures, with the MR 5 16 and 
MR 661 configurations, respectively, centered on the galaxy nu- 
cleus. Slight dithering between consecutive exposures was ap- 
plied to allow for a better rejection of bad pixels or faulty spec- 
tra. Details are provided in Table [1] 

All the OASIS data were reduced using the XOasis data re- 
duction softwareQ developed by CRAL (Arlette Pecontal and 
coll.), and adapted to the above-mentioned setups. This requires 
the creation of a spectral mask, following a detailed model of the 
optical path for the instrument, which allows an accurate optimal 
extraction of the signal from the raw exposures. The detailed re- 
duction pipeline includes CCD corrections (bias, fringing, re- 
moval of the overscans), wavelength calibration, flat-fielding 
(using both lamp calibrations and sky exposures), and flux cali- 
bration. Each individual exposure is then truncated to a common 
wavelength domain (from 4770A to 5530A for MR 5 16, and 
from 6620A to 6990A for MR 661). We then reconstruct images 



http://www-obs.univ-lyonl.fr/labo/oasis/download/ 
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Table 1. Characteristics of the SAURON and OASIS datasets. 



IFU 


Mode 


# Expo 


Exp. Time 


Spectral Sampling 


Spatial Sampling 


Expo FoV 


Total FoV 


OASIS 


MR 5 16 


5 


30 min 


2.15 A 


0.2" 


10.3" x7.4" 




OASIS 


MR 661 


4 


30 min 


2. 14 A 


0.2" 


10.3" x7.4" 




SAURON 


LR 


26 


30 min 


1.15 A 


0.8" 


42" x33" 


135" x90" 



by fully integrating over the spectral direction, and recenter each 
individual data cube so that the center of NGC 695 1 corresponds 
to (0, 0). The individual exposures are then merged, taking into 
account the effect of differential refraction, by re-projecting the 
cubes onto a common spatial grid centered at (0, 0) with spaxels 
of 0'.'2 X 0'.'2. The final merged cubes have a spectral sampling of 
2. 15 A for the MR 516 configuration, and 2. 14 A for the MR 661 . 

2.2. SAURON data 

NGC 695 1 was observed with SAURON at the WHT (iBacon et alJ 
12007 ) during a run in August 2003. Individual exposures cover 
a field of view of about 41" x 33", with spaxels of 0'.'94, a spec- 
tral resolution of 4.2 A and spectral sampling of 1 . 15 A over the 
range 4780 A to 5350 A . The spectral coverage of the SAURON 
datacubes overlaps with the OASIS MR 5 16, though it is signifi- 
cantly shorter in wavelength coverage. The coarser spatial sam- 
pling allows for a much larger spatial coverage of the disk of 
NGC 6951 than with OASIS, see also Table[T] We obtained a to- 
tal of twenty-six 30 min exposures, the goal being to cover as 
much of the stellar bar as possible. The total time spent on each 
pointing was adopted according to the average surface bright- 
ness, and included some slight dithering for each pointing to 
prevent systematics. The final area covered by the complete mo- 
saic is 135" X 90" (see Fig.[T]). The data fr om the single, centra l 
exposure has been presented previously bv lDumas et al.l(l2007h . 

The data reduction of all SAURON exposures was performed 
using a pipeline based on XSauron which is similar to the 
XOasis version, but adapted to the different optics and setup. 
The reduction steps are thus very similar to the ones described 
in the previous Section for OASIS. The biggest differences lie 
in the characteristics of the spectral mask, and the existence, for 
SAURON, of about 150 dedicated sky lenses which allow for an 
accurate sky subtraction for each exposure. 

The fully calibrated individual exposures were then recen- 
tered (using a direct image of the galaxy as a reference) and 
merged, by projecting each data cube on a common grid with 
O'.'S xO'.'S spaxels. To increase the S/N, the spaxels were Vonoroi 
binned to a S/N of 40. 

2.3. Spectral fitting of tlie IFU datasets 

We applied standard routines to extract the kinematics: the 
stellar kinematics (velocity, velocity dispersion and Gauss- 
Hermite and h4 val ues) with pPXF (£appellari & EmsellemI 
l2004t lEmsellem et alJ l2004) usin g a set of stell ar templa tes ex- 
tracted from the MILES lib rar>fl (ISanchez-Blazquez et all2006t 
iFalcon-Barroso et"ani201 ih and the IDL implementation avail- 
able from Michele Cappellari's web pag^E and the gas kinemat- 
ics via a version of GANDALF (kind ly adapted for our pur- 
pose by Marc Sarzi, ISarzi et aH l2006b . Fitted emission lines 
include U/3, the doublets [OIII] (A (4959,5007) A), [NI] (A 



^ http://miles.iac.es 

^ http://www-astro.physics.ox.ac.uk/~mxc/idl/ 



(5200,5202) A), in both SAURON, and OASIS MR 5 16, and [OI] 
{A (6300,6364) A), Ha {A 6562.8 A), and the doublets [Nil] 
(A (6548,6583) A) and [SII] (A (6717,6731) A) for the OASIS 
MR 661 data cube. The kinematics of the emission lines was 
fixed to the [OIII] doublet (SAURON) and [01], [Nil] and [SII] 
kinematics were fixed together in OASIS MR 661. The Ha emis- 
sion was strong enough to be fit independently. We emphasize 
that H/3 and Ha are present both in emission and absorption, 
however, the combined use of pPXF and GANDALF does take 
this into account. 

The resulting emission line maps from the spectral fitting for 
OASIS are shown in Fig. |2] It can be seen that the OASIS FoV 
just covers the circumnuclear ring. The Ha and H/3 emission is 
mostly concentrated in the ring, while the [OI] and [OIII] emis- 
sion is limited to the nucleus. The [Nil] and [SII] emitting gas 
is present in both the nucleus and ring. There are bright hot- 
spots in the ring that are seen in Ha, H/3, [Nil] and [SII]. The 
velocity field and dispersion are only shown for Ha, but are very 
similar for the other emission lines. The increase in the velocity 
dispersion in the western half of t he ring was also see n in the 
SAURON observations presented by iDumas et al.l (l2007l) and co- 
incides with the lowest intensities of Ha and HjS line emission in 
the ring. 

The spectral fitting results from the SAURON mosaic are pre- 
sented in Fig. [3] In the top-left panel of Fig. [3] we show the to- 
tal intensity over the FoV of the mosaic. The peak intensity of 
the central region has been saturated to show the details of the 
entire stellar bar region. The stellar bar is orientated approxi- 
mately east-west on the sky. The few high-intensity peaks are 
foreground stars. The stellar velocity field and dispersion are 
also shown in the top row. 

In the bottom row of Fig.|3]the H/3, [NI] and [OIII] emission 
line maps are presented. It is clear that their emission is con- 
centrated in the circumnuclear region. In the case of H/3, a keen 
eye can make out the ring shape of the emission in this figure 
(exactly as was seen in the OASIS map). 

2.4. HST images 

We retrieved F547M, F606W, F658N, F814W and F330W im- 
ages of NGC 6951 from the MAST HST archive. We extracted 
them from the archive calibrated with the "best" reference files. 
We also utiHze NICMOS FllOW and F160 W images from the 
HST archive al ready used in previous work (IVan der Laan et aP 
12011 . original: Hunt & Malkan (2004)). Finally, we obtained a 
Paa image, which was observed during the same campai gn and 
reduced using the same approach as for NGC 1097 in Hsi eh et al.l 
(1201 ll) . The full sample and detailed data reduction for the Paa 
image are presented in Liu et al. (in prep.). A list of all images 
used is given in Table |2] 

The Ha image was carefully continuum subtracted with 
the F606W image, as was the Paa (F190N) image using the 
F187N image. The F658N filter passband includes emission 
from the [Nil] double line emission. From spectra presented 
in IStorchi-Bergmann et all (l2007l) the ratio [NII]6583 A/Ha is 
known to be ~0.3 throughout the ring, implying a 37% contri- 
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Fig. 2. Top panels and bottom right: OASIS intensity maps of the emission lines Ha, Hj8, [OI], [OIII], [Nil] and [SII]. The color 
scale in each panel indicates the square of the intensity, the range of which is given below the color bar. This scaling was chosen to 
emphasize faint structures in the maps. The 'contact points', discussed in the text, are indicated with a star in the Wa panel. Bottom 
left panels: the velocity and velocity dispersion map of Ha. The flux contours of the integrated Ha emission are over plotted for 
emphasis. The velocity range in the bottom-left panel is [-150, 150] km s"', with an assumed systemic velocity of NGC 6951 of 
1425 km/s. The velocity dispersion range in the adjacent panel is [0, 200] km/s. The orientation of all maps is north up, east left. 
All positions are relative to the nucleus. 



Table 2. HST archival data 



Filter 


Program ID 


Exposure [s] 


Obs. Date 


PHOTFLAM 


Instrument 


Scale ["/pixel] 


F187N* 


11080 


1151.767 


13-02-2007 


3.6837500x10^'" 


NICMOS NIC3 


0.10 


F160W 


7331 


255.923 


16-12-1997 


2.3600094x10-" 


NICMOS NIC2 


0.076 


FUOW 


7331 


255.923 


16-12-1997 


4.3320170x10-" 


NICMOS NIC2 


0.076 


F547M 


5419 


300 


16-07-1994 


7.691 168x10""* 


WFPC2/PC 


0.046 


F606W 


8597 


160 


18-11-2000 


1.89971 1x10""* 


WFPC2/PC 


0.046 


F658N 


9788 


700 


24-10-2003 


1.9597158x10-'** 


ACS WFCl 


0.05 


F814W 


9788 


120 


24-10-2003 


6.9255574x10-20 


ACS WFCl 


0.05 


F330W 


9379 


1200 


17-08-2002 


2.2688470x10-'** 


ACS HRC-FIX 


0.025 



Notes: List of the observations extracted from the HST archive. *The Pao- image has been corrected for 
continuum flux with a F190N-filter image taken at the same time. 



bution to the observed flux from [NI1]6548 A and [NII]6583 A 
in the HST F658N filter for a V,,,, of 1425 km/s, for which we 
correct. 



3. NGC 6951 's circumnuclear ring so far... 

Before the start with the analysis of the stellar properties of 
NGC 695 1's circumnuclear ring, a brief summary of the gaseous 
properties of the circumnuclear region that are known to date is 
presented. For clarity later on, all numbers pertinent to the dis- 
cussion later are also listed in Table |3] The large scale stellar bar 
in NGC 6951 has a PA of 84° on sky (Mulchaey & Regan 1997). 
Its semi-major axis is 3.0kpc. The dust lanes of the bar are thus 
orientated almost horizontally on the sky. The contact points in 
the ring are in the north and south of the ring, see also Figs. |2] 
andg] 

'^CO(l-O) and '^CO(2-l) observations, taken with the 
IRAM PdBl and 30m telescop e, of the circumnuclear region 
of NGC 6951 are presented by iVan der Laan et al.l (120111) . The 



resolution of the CO(l-O) observations is 3.11" by 2.59", the 
CO(2-l) observations are 1.72" by 1.56". The FoV of the CO(l- 
0) observations is shown in Fig.lT] the FoV of the CO(2-l) data 
is about half that size. Both completely cover the ring. Their 
r.m.s. noise corresponds t o M//2 ~ 10^ Mo/beam. We refer to 
I Van der Laan et alT(l201 ih for further details on the data reduc- 
tion. A total H2 molecular gas mass of 2.2 x 10'' M© was derived 
from the CO(l-O) observations, assuming a Galactic conversion 
factor and corrected for helium abundance. An H CN(l-O) in- 
tensity map obtained with PdBl was presented by iKrips et al] 
(2007). A map with slightly different Mv-weighting from that 
dataset was kindly provided by M. Krips (priv. comm.), with a 
beam size of 1.89"by 1.83"(PA 78.7°). 

The key points for the molecular gas morphology are promi- 
nent gas spiral arms, that are co-spatial with the dust lanes along 
the large-scale stellar bar and end at the radius of the circumnu- 
clear ring. The CO(2-l) observations, which have a higher spa- 
tial resolution, show that the spiral arms continue into a ring. The 
HCN(l-O) traces higher density molecular gas than CO, and the 
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Fig. 3. SAURON intensity maps of the stellar light (top left), H/3 (bottom left), [NT] (bottom middle) and [OIII] (bottom right) emission 
lines, as well as the stellar velocity field (top middle) and velocity dispersion map (top right). The color scale in each panel indicates 
the square of the intensity, or the velocity, the range of which is given below the color bar The orientation of the maps is north up, 
east left. All positions are relative to the nucleus. To guide the eye, the contour corresponding to the central H/3 emission has been 
over-plotted onto the velocity and velocity dispersion plots. The velocity range in the top-middle panel is [-200, 200] km s the 
assumed systemic velocity of NGC 6951 is 1425 km/s. The velocity dispersion range in the top-right panel is [0, 185] km/s. 



Table 3. Properties of NGC 695 1 



Parameter 


Value 


Reference 


Type 


SAB(rs)bc 


(1) 


Nuclear Activity 


LINER/Seyfert 2 


(2) 


Dynamical Center (Locus Radio Continuum) 




RA (J2000) 


20''37"'14.123' 


(3,4) 


Dec (J2000) 


66° 06' 20.09" 


(3,4) 


Inclination Angle 


46.2° 


(5) 


Position Angle 


138° 


(5) 


Adopted Distance 


24.1 Mpc 


(6) 


Stellar bar PA 


84° 


(7) 


Mass mol. gas (r<3 kpc) 


1.6 xlO' Mo 


(8) 


Mol. gas inflow rate onto rinj 


; 2.0Mo/yr 


(8) 


Large scale bar pattern speed 


28.9±1.2km/s/kpc 


(9) 


V 

' sys 


1425 km/s 


(5) 


Period pattern speed 


210Myr 


(9) 


1 rotation for a ring orbit 


24Myr 


(9) 



Reference s: (1) [ de Vaucouleu rs et al 
(3) LE DA, iKripset al. ( 2007), '~ 
i)|Tullv I 



li2009l). (6)|Tu11v (1988), (7) Mu 
IVan der Laan eral.rj20Tll) . (9) this work 



19911). ( 2) iPerez et al.l ll2000h 

(4) Saikia et al. (2002). (5) Ha an et alj 
Mul chaev & Regan. (.199% . (8) 



HCN intensity distribution shows a similar moiphology to the 
CO(2-l) molecular gas. However, th e intensity maxima are at 
slightly different positions in the ring. IVan der Laan et all (1201 ll) 
derived, from modeling the orbits in the circumnuclear region, 
that the inflow rate of gas onto the ring is ~2.0 M© per year. 

The HI distribution, presented in Haa n et al.l (l2009l) . shows 
that the circumnuclear region of NGC 695 1 is Hl-deficient. The 
neutral ISM in the circumnuclear region is thus molecular gas 
dominated. Given an estimate for the dynamical mass at 1 kpc, 
(Mdyn = ^ ^ 8 X 10'^ Mq, with v=200km/s), atomic and molec- 



ular gas makes up approximately 25% of the dynamical mass 
within this radius. 

Two determinations of the large scale bar pa ttern speed of 
the stellar bar in NGC 695 1 exist in the literature. lAguerri et al.l 
(1998) obtained a value of Q.p - 30.3 km/s/kpc (when compen- 
sating for the difference in adopted distance) based on the rota- 
tion cu rve and location of the corotation radius and iHaan et aP 
(I2009h derived a value of Qp = 23 km/s/kpc based on an anal- 
ysis of dynamic resonances. The SAURON IFU mosaic cov- 
ers most of the large scale stellar bar of NGC 6951, and is 
ideally suited to measure the pa ttern speed using the model - 
free Tremaine- Weinberg method (iTremaine & Weinberglll984l) . 
Using the SAURON mosaic gives a value of 28.9+1.2 km/s/kpc, 
coiTected for inclination (46.2°) and at a distance of 24.1 Mpc. 
Details will be available in a forthcoming paper 

With a pattern speed of 28.9 km/s/kpc, the pattern will take 
210 Myr to complete one rotation. The gas and stars rotate with a 
velocity of ~150km/s at the ring radius. They will thus complete 
a rotation in the ring in 24 Myr 



4. Extinction in the ring 

From the three Hydrogen recombination lines (Ha, H/3 and Paa) 
in our data sets we can compute a spatial distribution of the ex- 
tinction in the ring. While determining the extinction is not di- 
rectly crucial to the overarching goal in this work, it provides a 
consistency check on the deiived color excesses of the star clus- 
ter fitting (Sect. |5]l. We use two different ratios to compute the 
extinction, the Ha/H/S ratio (as derived from our OASIS data) and 
the Ha/Paa ratio (as derived from the HST data). 

Both the Ha and li/3 maps are obtained during the spec- 
tral fitting of the OASIS data cubes. The fitting takes Ha and 
H/3 stellar absorption into account as explained in Sect. 12.31 In 
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The Ha/HjS based color excess map (OASIS) and the HafPaa 
based color excess map (HST) show similar distributions, and 
have similar values in the range [0.2 - 0.5] mag. 



5. Star clusters 




1QD pc 



Fig. 4. False color image of the circumnuclear ring in NGC 695 1 , 
based on the F814W (red), F606W (green) and F547M (blue) 
broad band filters. The image (~ 14" x ~ 14") is orientated 
with north pointing up and east to the left. The center has been 
saturated for better visibility of the ring, and bad pixels have 
been removed. The 'contact points' of the ring are indicated with 
blue circles. 



the c ase of case B recombination, the ratio Hc/HyS will be 2.86 
jOsterbrock 1989). Extinction will increase the value. In indi- 
vidual pixels, we find ratios as high as ~4, which is equal to an 
extinction of ~ 1 magnitude. Inside the ring, the surface bright- 
ness of especially Hfi is so low, that no reliable Hafli/3 ratio can 
be computed. Where that is the case, the position has been left 
blank. In Fig. |5] (right) the color excess, E(B-V), based on the 
Ha/HjS ratio is plotted. Color excess is a wavelength independent 
measure of extinction, and is related to extinction via E(B-V) - 
Ay/Ry, where Ry gives the total-to-sele ctive extinction, wh ich 
we assume to be 4.05 (starburst value of lCalzetti et al.llToOOh . 

The FoV of OASIS only barely covers the entire ring. 
Nevertheless, it is possible to see some trends in the Hafli/3 map. 
The extinction is clearly higher near and 'down-stream' from the 
contact points between ring and bar-induced gas lanes (galaxy 
rotation is clockwise). 

The ratio of the HST Ha and Paff observations provides an- 
other measure of the extinction. The advantages of these ob- 
servations are the complete view of the ring and higher spatial 
resolution (Fig. |5] lef t). Without extinction the expected ratio 
HafPaa is about 8.5 (iHummer & Storevlll987l case B recom- 
bination). Inside the ring the intensity of both Ha and Paa is 
too low to obtain reliable ratios. Those positions have been left 
blank. Again we see the highest extinctions 'down-stream' from 
the contact points. In the western part of the ring we also find 
'pockets' of higher extinction which can be related to the highest 
intensity CO and HCN peaks. These places are heavily obscured 
and dense, and thus are likely the sites of future star formation. 
Comparison with the Ho' distribution shows the brightest Ha re- 
gions to be anti-correlated with extinction. 



Star clusters are a prime way to study the ages of stellar pop- 
ulations; all the stars within a cluster are coeval and clusters 
are bright and compact, thus easily observable even over large 
distances. At the distance of NGC 6951 star clusters are point 
sources and can easily be separated from the background light 
in the galaxy. 

5.1. Identification 

Two independent routines were used for the detection and pho- 
tometry of the star clusters in the ring of NGC 695 1 . These rou- 
tines are the DAOPHOT package (Stetson '1987"; Stetso n et al. 
1990; Stetson 1992) in IRAF and SExtractor (Bertin & Amouts 
1996h . With this approach, the strengths of both methods can be 
exploited to obtain the most complete star cluster candidate list. 

For the detection of the star clusters a composite optical 
image was created out of 4 individual optical bands; F547M, 
F606W, F814W and F330W. First, the F606W, F814W and 
F330W images were aligned with the F547M image using stan- 
dard procedures in IRAF For the F814W and F330W bands the 
images were resampled to the F547M pixel size while conserv- 
ing the flux. The four resulting images were then normahzed and 
added together 

Applying DAOFIND, 33 star clusters candidates were de- 
tected in the composite optical image (within r ~20"). The star 
cluster detection in the NIR images was done separately from 
the composite optical image, as the NIR images have a lower 
angular resolution. The FllOW image was used to detect the 
NIR star clusters and the resulting positions were applied to the 
F160W image. Most of the NIR detections have a counterpart in 
the composite optical image, but several sources were detected 
without one. This leads to a total of 42 detected star clusters can- 
didates using the DAOPHOT method. 

Next, SExtractor was run. The strength of SExtractor is its 
ability to 'de-blend' crowded regions. As can be seen from the 
false color image (Fig. HJ, especially the northern region of the 
ring is very crowded. For detection and the subsequent photom- 
etry of star cluster locations the same composite optical image 
from the DAOFIND analysis was used. SExtractor automatically 
determines apertures around the selected source. If the individ- 
ual images were used independently, the selected apertures be- 
tween the different wavelengths could vary, which would intro- 
duce systematic changes in our photometry. By using the com- 
posite optical image as the 'detection image' this is avoided and 
each source will have the same aperture at the different wave- 
lengths. 

The SExtractor output is filtered by discarding all detections 
flagged as 'bad' by SExtractor itself (e.g. likely to be cosmic rays 
or one or more pixels that are saturated). The likely candidates 
are further constrained by location (a rectangular area around 
the ring of 12" x 12" ) and by relative error in the measured 
flux (<20%). 43 star cluster candidates are detected. 

In total, 55 different star clusters candidates are detected (see 
Table Ell, about half of which are detected by both DAOPHOT 
and SExtractor (28 out of 55 or 5 1 %). 
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Fig. 5. Left: Color excess E(B-V) in the ring as measured from the Hq,/H^ line ratio using the OASIS line maps. Right: Color excess 
in the ring as measured from the HafPaa line ratio using HST images. The pixel scales of each map are given in the lower-right 
comer of each panel. 
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Fig. 6. Locations of the candidate star clusters as detected with 
IRAF/DAOPHOT and SExtractor (open circles). W e have also 
indicated the locations of the star clusters found by iBarth et aP 
(1995) (grey squares). The axes denote the projected distances 
in arcseconds from the photometric nucleus, see Table |3] 



5.2. Properties 

The 55 st ar cluster candidates can be compared with the work by 
iBarth et al. (1995), who used a single F547M image (the same 
one we extracted from the archive). All 24 sources already pre- 
sented in that work are also included in our list. The additional 
candidates lie in the northern region and several have lower sur- 
face brightnesses. Investigation of the F330W image, which has 
the finest native pixel scale, shows that some of our detections, 
particularly in the crowded northern region of the ring, are actu- 
ally still composed of two or multiple point sources/star clusters. 



SExtractor simultaneously detects the star cluster candidates 
and measures their flux. For the star cluster candidates identi- 
fied with DAOFIND, their photometric fluxes were obtained us- 
ing PHOT (also in the IRAF/DAOPHOT package). For the mea- 
surement the F330W and F814W images are used at the reso- 
lution of the F547M and F606W images. In PHOT die circular 
apertures for the optical images were chosen to have a radius 
of 2 pixels (0.09")- This radius contains most of the observed 
PSF (4 pixels x FWHM ~ 20 pc), while minimizing the contam- 
ination from other sources (such as other close-by star clusters). 
The sky subtraction was based on an annulus around each source 
with radii between 6 and 16 pixels. Flux outliers in this annulus 
(e.g. other clusters or cosmic -rays, which would over-estimate 
the sky-levels) were automatically excluded from the sky esti- 
mate. For the NIR images an photometry aperture radius of 1.2 
pixels was used, which corresponds to a size equal to the optical 
aperture. The radii of the sky estimate annulus were again 6 and 
16 pixels for the NIR images. 

The Tiny Tim packagqj was used to obtain the necessary 
aperture corrections for each band. The aperture corrections are 
for the F330W band 1 . 16 x/330, F547M band 1 .29 x fs^, F606W 
band 1.31 X/6O6, F658N band 1 .40x/658, F8 14W band 1 .36x/8i4, 
Fl low band 1.38 x /no and F160W band 1.59x/i6o. 

All measured photometric fluxes are converted i nto absolute 
magni tudes, using a distance modulus of 31.9 mag. lBarth et al.l 
(Il995h used a distance modulus of 3 1 .2, which corresponds to 
a distance of 17.4 Mp c. Assuming a dist a nce o f 24.1 Mpc to 
NGC 6951, as done in IVan der Laan et al.l (1201 Ih increases the 
distance modulus to 3 1 .9, or correspondingly increases the mag- 
nitudes by -0.7 mag. The absolute magnitudes are given in Table 

a 

The absolute magnitudes from b oth methods are in good 
agreement with each other and with Barth et al ] (ll995h . with at 
most a 0.5 mag difference in the optical bands and I mag in the 
NIR bands, mostly due to differences in sky subtraction. These 
values have not (yet) been corrected for extinction. A measure 
for the intrinsic extinction is derived during the star cluster fit- 
ting. 



http://www.stsci.edu/software/tinytim/tinytim.html 
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Fig. 7. Results of the;i'^ fitting. The four segments of the ring are indicated in the plots. From left to right are shown: stellar age, 
mass and color excess of the cluster The individual results are averaged every 45°(crosses) and the variance is indicated with the 
error bars. The angle starts at the southern point of the ring and is measured counter-clockwise. 



6. Ages of the stellar populations in the ring 

Determination the ages of the stars in the ring gives information 
on the age of the ring in NGC 695 1 itself, since the star formation 
occurs in the circumnuclear ring. 

6.1. Starclusters 

The onHne Starburst99 (v6.0.2) code ("Leitherer etaD 119991; 
["Vazquez & Leitherer 2005; Leitherer et al. 2010) was run to ob- 
tain model SEDs for star clusters of different ages. We assumed 
a fixed stellar mass of 10^' Mq, an instantaneous star formation 
event with a Kroupa IMF and Padova AGB solar metallicity 
stellar tracks. The model SEDs span the UV to the NIR in the 
time range of lO'' to 10^ years in steps of 2x10'' years. We also 
include nebular emission in the model fluxes, which adds sig- 
nificant emission for young ages at, predominantly, the shorter 
wavelengths. 

For comparison to our observed data points, these model 
spectra are weighted with a range of E(B-V) values ([0, 3] in 
0.01 steps) and stellar masses ([10^,10^ Mq], in 0.2 dex steps). 
Mass scales linearly with luminosity and as such scales to the 
observed SED, independent of wavelength. A different distance 
modulus (31.2 instead of 31.9) would have a similar effect. To 
convert the c olor excess int o an extinction curve the Calzetti ex- 
tinction law (ICalzetti et al.l l2000) for starburst galaxies is used. 
The Starburst99 spectra so obtained are convolved with the dif- 
ferent HST filter throughput curves, and converted to magni- 
tudes. A x^-fit is applied to each star cluster set of observed 
magnitudes. The lowest value is selected as the best fit to 
our data (Table |5]l. 

Table |5] shows that the best fit ages for the star clusters span 
young (4 Myr) to intermediate (300 Myr) ages, with several older 
clusters present (^700 Myr). The stellar mass range of the star 
clusters is ~10'*" to ~10''-^ Mq and the E(B-V) color excess 
varies from to 1.0 magnitudes, with most star clusters hav- 
ing a color excess of <0.6. The results of the star cluster fitting 
are also shown as a function of azimuthal position in the ring in 
Figure 12] The four quadrants (SW, NW, NE and SE) should be 
visualized as overlaying compass points on the ring and are cho- 
sen with the contact points between gas/dust lanes of the large 
scale bar and the ring in mind. These points fall just inside the 
SW and NE quadrants. The rotation sense of the galaxy is clock- 



wise, which corresponds to an decreasing angle. The results are 
averaged every 45° to detect any azimuthal trends in age, mass 
or color excess. 

The star clusters have predominantly intermediate age and 
are massive, which explains why they are still detectable. From 
the left panel (age) in Fig. [T] it can be seen that the star clus- 
ter ages in the ring vary little, although they do drop in the SW 
and NE quadrant, before risin g. It is not clear wheth er this trend 
could be a weak age gradient. Mazzuc a et al. I (l2008h find a qual- 
itatively similar age distribution, although we should point out 
that ~ 100 Myr ages for the Ha emission region in that work 
can only be achieved with a continuous star formation history. 
The star clusters are assumed to have for med in a singl e burst . 
However, the Ha regions discussed in Mazz uca et al ] (I2008h 
cover a much larger area than the individual star clusters in this 
work. Thus, both results would be in agreement if some younger 
stellar population is not or no longer concentrated in star clus- 
ters, but is present in the ring. We discuss this scenario further in 
the next subsection. 

The color excess derived from the star cluster fitting (on aver- 
age 04-0.6 mag) is somewhat larger than the color excess values 
obtained from the Ha/Hfi and Ha/Paa ratios presented in Sect. |4] 
(on average 0.4 mag). However, the hydrogen emission and the 
star clusters do not come from the same physical locations. The 
extinction that is found for the star clusters is confined to the 
star cluster region, and can thus be higher than the more global 
extinction map derived from the Hydrogen lines. 

6.2. Stellar populations from SAURON IFU 

The star clusters discussed above (Sect. 6.1) have ages of sev- 
eral hundred Myr and are visible due to their concentrated light. 
How do the star clusters relate to the underlying diffuse stellar 
populations in the ring? Does the ring have a different star forma- 
tion history from the disk, or does the disk contain similar stellar 
populations? With the SAURON IFU dataset these questions can 

be addressed. 

Similar to the work bv lStoklasova et aP (l2009 h. we fit a com- 
bination of 'young' and 'old' single stellar population (SSP) 
spectra to each Vonoroi cel0 of the SAURON mosaic using pPXF 



' At the ring positions the size of a Vonoroi cell is equal to the origin 
spaxel. 
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Fig. 8. The stellar 'mass' maps in [a.u.] of mass by age bin obtained by fitting the SAURON IFU dataset with a combination of 4 SSP 
spectra. The maps are corrected for inclination and are oriented north up, east left. The approximate location (a- 5", b~ 4.5") of 
the ring is indicated by the dashed ellipse. The 'mass' uncertainty is within 10%. 



and GANDALF with similar settings as in Sect. 12.31 The goal is 
to split the light at each position into its 'young' and 'old' stel- 
lar population fractions. The assumption here is that the stellar 
content observed in each cell can be approximated by a finite, 
small, combination of SSPs. Due to the continuous evolution of 
stars over time, this assumption can usually be made, as long as 
the combination is chosen to remain sensitive to the full possible 
age range (0 - 13Gyr). SSP spectra sp anning - 13Gyr were 
extracted from the MILES datab ase^ jSanchez-Blazquez et alj 
120061: iFalcon -Barroso et al ] l2011h with a spectral resolution of 
cr - 60km/s and spectral sampling of 1. ISA/pixel. In order to 
avoid an age-metallicity degeneracy, only solar metallicity SSP 
models were selected. The original fit of the SAURON data, dis- 
cussed in Section |23] was more extensive; a range of metallici- 
ties was used there. It showed that solar metallicity is the dom- 
inant metallicity and thus a valid simplification here. Then, us- 
ing mock observations, generated from the MILES SSP spec- 



* http://miles.iac.es 



tra together with added noise, the set of SSP spectra that still 
reasonably span the full possible age range could be restricted 
to 4 SSP spectra. Their age identifiers in the MILES database 
are 70.8Myr, 158.5Myr, 398.1Myr and 3.2Gyr. From here on, 
they will be designated 'young', 'intermediate I', 'intermediate 
ir and 'old'. These 4 spectra will obviously not give the same 
age refinement as the full set, but our main goal, as stated, is 
to separate 'young' and 'old' stellar light at each spaxel/Vonoroi 
cell. Each Vonoroi cell in the SAURON mosaic was fit with a com- 
bination of these 4 SSPs. 

The spectral fitting routine pPXF gives a (normalized) num- 
ber weighted output of the SSPs used to fit the spectra at each 
position. We used a 'bootstrap' -like method to obtain error es- 
timates on the output generated by fitting the 4 SSPs to each 
spaxel/Vonoroi cell. From the fit of the data with a full set of 
SSP models, a measure for the fitting residual was obtained. 
This 'residual', which is approximately 3% of the flux, was used 
to introduce random noise of this magnitude into the spectra. 
The fitting was iterated 30 times, each time introducing new 
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random noise. Again based on mock data, the uncertainty of 
the number weighted fractional contribution by each of the 4 
SSPs is estimated to be within 10% after 30 iterations. We no- 
tice from the mock data that if a large 'old' population is present, 
the fitting underestimates the 'young' fraction by several per- 
cent. The opposite is not true. The results therefore remain sen- 
sitive to older populations also in regions dominated by light 
from 'young' stellar populations. Further, the intermediate frac- 
tions are slightly under-predicted in all cases, to the advantage 
of the 'young' and 'old' fractions. The fitting routine attempts 
to minimize the number of SSPs necessary where possible. The 
normalized number weighted output is transformed into (again 
normalized) luminosity fractions, with the known luminosity of 
each SSP spectrum. The luminosity fractions can also be trans- 
formed into normalized stellar mass fractions, by use of the M/L 
ratio tabulated for each SSP spectrum. Finally, we multiply the 
normalized mass fractions with the stellar intensity at each cell 
(Fig. [3] top-left panel) to obtain a common baseline for all spax- 
els. This can be interpreted as an 'absolute' mass as function of 
spatial position for each age bin. The resulting maps for the 4 
age bins are shown in Fig. [8] 

As can be seen from Fig. [8] (results have been corrected for 
inclination), the ring is prominent in the 'young' bin, compris- 
ing about 30-40% of the mass at its galactocentric distance. The 
contact-points between the large scale bar induced gas lanes and 
the ring stand out at the north and south positions in the ring. 
The mass of 'young' stars is highest at those positions. In the 
'old' bin the opposite picture is observed; almost no 'old' stel- 
lar population is clearly associated with the ring, the 'old' stars 
are distributed in a typical bulge-like profile. Significant mass in 
'intermediate' age stellar populations is found at the inner edge 
of the ring. 

Fig. [8] shows that the ring is not circular, but slightly el- 
lipsoidal, when corrected for inclination. Ellipticity is very 
common for the entire pop ulation of circum nuclear rings (e.g. 
iButa & Crocker.. 1993: Comeron et al.ll2010l) . 

7. Discussion 

7.1. Star formation signatures in fhe ring 

By combining the results from the star cluster analysis and the 
underlying stellar 'field' population through the IFU fitting, the 
stellar content in and around the circumnuclear ring has been 
fully constrainted. What does it say about the manner in which 
star formation proceeds in this circumnuclear ring? Is there an 
azimuthal ag e gradient, as in the 'pearls-on-a-string' scenario 
proposed by Boker et al.l (l2008h . or has star formatio n resulted in 
a hom ogeneous ('popcorn') stellar age distribution (tEbnegreenl 
[T99I ? 

The recent star formation episodes are concentrated in the 
ring region (see Fig. |8]l. This is in accordance with our open- 
ing statement that the circumnuclear ring is the gas reservoir for 
the circumnuclear region. The 'old' stellar bin in Fig. |8]mainly 
shows a bulge-like profile, with no additional mass at the ring 
radius itself. The IFU fitting shows a significant 'young' mass 
fraction in the ring itself, between 30% and 40%. The aver- 
age age of the star clusters in the ring is young/intermediate; 
200-300 Myr This discrepancy in age (range) between the two 
methods is partly due to the IFU method, which slightly under- 
predicts the intermediate stellar fraction in favor of the 'young' 
fraction. However, the IFU fitting also covers the full ring and is 
thus equally sensitive to star formation which did not lead to the 
formation of > IQ'^Mq star clusters (the lightest clusters in our 



sample). The ring could be currently in a state where star forma- 
tion leads to < 10"* M0 clusters. Possibly, the current gas inflow 
rate is lower than > 100 Myr ago when the detected star clusters 
were formed. 

Both the star cluster and IFU results show that the distri- 
bution of the 'young' stars is not uniform throughout the ring. 
The relation between the contact points and star formation is the 
main observable for the 'pearls-on-a-string' scenario for star for- 
mation. The locations with the highest mass in 'young' stars are 
close to the gas inflow points, the contact points, which is a point 
in favor of the 'pearls-on-a-string' scenario. 

In Fig. |9] we further plot the intensity distributions of the 
molecular gas tracers CO(2-l) and HCN(l-O) as a function of 
azimuthal angle in the ring. The azimuthal distribution of the 
'young' stellar mass fraction (in percentage of the total mass 
at each position, based on the IFU results) is also plotted there. 
The highest concentrations of molecular gas are close to the con- 
tact points in the north ((p ~ 180°) and to a lesser extent in the 
south (0 ~ 350°). The variations are especially pronounced in 
the HCN distribution. The fraction of the total mass in 'young' 
stars by comparison shows much less variation. 

The highest concentrations of 'young' stars are found in the 
NW and SE segments of this figure. The wave pattern of the 
large scale stellar bar is slower (28.9km/s/kpc pattern speed) 
than the stellar and gas velocity (~ 150km/s) at the circumnu- 
clear ring radius (580 pc). This means that the overdensity points 
will move counter-clockwise in the inertial frame of the galaxy 
(galaxy rotation is clockwise). If stars form preferentially at or 
near the overdensity points, the stars and gas will move faster 
and out of the overdensity region. The fraction of 'young' stellar 
mass must thus be higher in front of the overdensity region (in 
Fig. |9]'in front of equals at lower azimuthal angle). The off'- 
set between the highest fractions of 'young' stars and the peak 
in HCN emission in the north is ~ 20°. Approximately 1.5 Myr 
would be needed to travel 20° at the ring radius. 

A half (two-fold symmetry) orbit of the bar wave pattern in 
the ring takes approximately 100 Myr, and stars that have formed 
at those locations will make 5 full orbits around the ring before 
the pattern rotation is complete. Hence, any age signature will 
be severely smeared within that 100 Myr pattern time, and erase 
the 'pearls-on-a-string' scenario in favor of the 'popcorn' one. 
The small indications of young fractions we find are thus all we 
can detect. 

7.2. Age of the circumnuciear ring 

The stellar age content of the ring can assist in the determination 
of the age of the ring as a dynamical structure. Taking together 
the stellar ages, gas inflow rate and dynamical times allows both 
a minimum and maximum age to be constrained. 

The lack of extra mass in the ring in the old stellar bin in Fig. 
[8]suggests that the ring cannot be older than ~3 Gyr The average 
age of the star clusters is 2 00-300 Myr with some older (up to 
1 Gyr). 1^ der Laan et al.l ([201 1.) already inferred a minimum 
age of the large scale stellar bar and thus the ring of 1 . 1 Gyr, 
based on mass inflow and total molecular gas mass (after heUum 
correction). 

The inflow rate determined in i Van der Laan et al 1 (1201 lb was 
2.0Mo/yr The total mass in the star clusters we found is ~ 
6xlO^M0. A star (cluster) formation rate history, based on the 
star cluster results, is given in Fig. [10] The highest SFR was 
reached 800 Myr ago, after which the ring was quiescent for 
300 Myr Star formation has started again in the last 400 Myr 
Over 1 Gyr (the oldest measured star cluster age), the average 
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Fig. 9. Azimuthal mass fraction of young stars (asterisk), combined with the azimuthal intensity distribution of CO(2-l) (red Hne, 
dashed +1 standard deviation) and HCN(l-O) emission (blue line, dashed +1 standard deviation). The angle corresponds to the 
south, angle increases counter-clockwise. 



star (cluster) formation rate was ~ O.O6M0/yr A star formation 
efficiency of a few percent is entirely plausible. That implies 
that the current inflow rate is a good estimate of the average 
inflow rate over the last 1 Gyr, even though the star formation 
rate was variable, and that the previously inferred minimum age 
of 1 . 1 Gyr, based on that inflow rate, is still plausible. If the star 
formation efficiency is higher, a lower inflow rate would suffice 
and the minimum age would be higher. 

The fraction of angular mom entum lost by gas per orbit, 
dL/L, is 0.2-0.4 (iHaanet al.1 12009) . which means gas needs at 
least 3-5 orbits to settle on the ring, which equals ~200 Myr (or- 
bit times are larger at larger radii). 

The oldest likely age for the ring is 3 Gyr, given the lack 
of a significant number of 'old' age stars from the analysis in 
Sect. 16.21 If the ring were 3 Gyr, that implies that the ring was 
quietly accumulating gas for up to 2 Gyr before star formation 
in the ring commenced. Assuming a constant mass inflow of 
2.OM0/yr and a ring of 600pc (5") radius and 100 pc (1") width, 
that would imply an absurdly high molecular gas surface density 
of 11 ,000 Mo/pc^ at the onset of star formation. That is much 
higher than the average lOO Mp/ pc^ value of most giant molec- 
ular clouds (e.g. Bigiel et al. 2008). To reach an average surface 
density of 100 Mo/pc^ would thus only take 19 Myr. 

The current age of the circumnuclear ring is therefore most 
likely between 1 .0 Gyr and 1 .5 Gyr 
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Fig. 10. Star formation rate history in the circumnuclear ring. 
The masses of the star clusters are grouped in 100 Myr age bins 
and divided by the width of the age bin. No correction has been 
applied for cluster disruption, which would increase the SFR at 
older ages. 



8. Conclusions 

In this work the properties of the stellar populations of the cir- 
cumnuclear ring in NGC 6951 were investigated. The goals were 
to estimate the age of the circumnuclear ring itself, to determine 
the age of the stellar population(s) present in the circumnuclear 
ring, and to determine whether the 'popcorn' or 'pearls-on-a- 
string' scenario holds true for the nature of star formation in the 



circumnuclear ring. The main results can be summarized as fol- 
lows: 

The ring is most likely between 1 .0 Gyr and 1 .5 Gyr old. It 
has a stellar population distinct from the rest of the circumnu- 
clear region that was formed over that time frame. The maximum 
possible age for the ring is 3 Gyr However, it is very unlikely 
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that the ring has been in place for that long without forming stars 
for its first l.SGyrs. 

Analysis of the general stellar content of the ring, by means 
of the SAURON mosaic, shows that 30-40% of the stellar mass in 
the ring is from 'young' (<100Myr) stars. The largest fractions 
are found near the contact points that link the ring with the large- 
scale stellar bar induced gas lanes. 

The majority of the star clusters in the circumnuclear ring 
have intermediate ages. The youngest star clusters (~100- 
200 Myr), when the results are averaged over 45 ° intervals, are 
in the north and south, close to the contact points were gas flows 
onto the ring. Such a connection is indicative of the 'pearls- 
on-a-string' scenario for star formation in circumnuclear rings. 
However, the dynamical time of the ring is short; it takes only 
24 Myr for a star cluster to complete one rotation of the ring, 
while the age gradient can only be renewed every 100 Myr. 
Therefore we do not find definite evidence of the 'pearls-on-a- 
string' scenario in the circumnuclear ring of NGC 695 1 . 
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Table 4. Star clusters identified with IRAF/DAOPHOT and SExtractor 
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1 "2 "2 

-ID.D 


A A1 
U.UI 


n A 
- /.U 


A 1 

U.21 


T A 

- /.U 


'2T 

2.D / 


C A 

-D.4 


A AT 

u.y / 


-J.J 


A '2 1 

o.Dl 


A n 

-4.y 


0'2 

2. 2D 




1 

1.0 


1.8 


n A 

-/.u 


r\ A 
U.48 


-ID. / 


A A 1 
U.UI 


n 

- 1 .1 


1 '30 
1.D2 


C 
-8.D 


1 AC 

1.46 














o c 

Z3 


o c 
2.D 


1.0 


1 A 
-iU.J 


A A'3 

U.UD 


1 A /I 

-1U.4 


A A/; 
U.Uo 


1 A '3 

-lU.D 


A A/1 

U.U4 


1 A C 
-lU.D 


A /I T 

u.4/ 


1 A A 

-IU.U 


All 
U.ll 


A A 

-y.y 


A '3 A 

U.DU 


A 

-y.2 


A '3A 
U.DU 




9 (S 
z.. u 


1 f, 
i .u 


Q Q 
-y.y 


OzL 




08 
u.Wo 


-10.1 


0.07 


-10.4 


w.uw 


1 


W. It- 


Q 7 


AVt 

W.T-W 


Q 4 


yj.jo 


26 


2.5 


1.1 


-9.4 


0.06 


-9.9 


0.07 


-9.8 


0.05 


-10.2 


0.56 


-9.7 


0.13 


-9.2 


0.41 


-8.5 


0.42 


27 


-3.6 


1.0 


-8.8 


10.00 


-9.4 


0.11 


-9.4 


0.09 


-10.6 


0.55 


-8.9 


0.23 










28 


-3.8 


0.5 


-9.9 


0.04 


-10.4 


0.06 


-10.3 


0.04 


-10.9 


0.40 


-10.2 


0.10 


-9.2 


0.41 


-8.5 


0.42 




-3.8 


0.5 


-9.7 


0.05 


-10.5 


0.07 


-10.5 


0.05 


-11.0 


0.46 


-10.3 


0.12 
























continued on next page 

















Table 4 - continued from previous page 



ID 


Aa ["] 


A6 ["] 


M330 


O"330 


M547 


0-547 


M6O6 


0-606 


M658 


0-65^ 


M814 


O" 814 


Mi 10 


o-uo 


Mi6() 


O" 160 


29 


-3.7 


0.3 


-11.1 


0.02 


-10.9 


0.05 


-10.7 


0.03 


-10.8 


0.42 


-10.4 


0.10 


-8.8 


0.49 


-7.6 


0.62 




-3.7 


0.3 


-11.2 


0.02 


-11.1 


0.05 


-11.0 


0.04 


-11.3 


0.39 


-10.6 


0.10 


-9.6 


0.40 


-9.0 


0.45 


30 


3.2 


-0.3 


-10.7 


0.02 


-10.6 


0.05 


-10.7 


0.03 


-12.4 


0.21 


-10.1 


0.11 


-8.9 


0.48 


-8.0 


0.54 




3.2 


-0.3 


-10.2 


0.03 


-10.1 


0.08 


-10.2 


0.06 


-11.9 


0.29 


-9.8 


0.15 


-9.1 


0.52 


-8.3 


0.68 


31 


0.2 


-0.3 


-9.6 


0.06 


-10.9 


0.05 


-11.5 


0.02 


-13.2 


0.14 


-11.3 


0.06 


-11.4 


0.15 


-11.2 


0.12 


32 


-5.5 


-0.4 


-9.3 


0.04 


-10.4 


0.06 


-10.4 


0.04 


-10.4 


0.50 


-10.2 


0.10 


-8.9 


0.46 


-7.9 


0.55 




-5.5 


-0.4 


-9.5 


0.04 


-10.7 


0.06 


-10.7 


0.04 


-10.8 


0.50 


-10.5 


0.10 


-9.4 


0.43 


-8.7 


0.48 


33 


-4.9 


-0.5 


-9.2 


0.06 


-10.2 


0.08 


-10.2 


0.05 


-10.2 


0.66 


-9.9 


0.14 


-9.0 


0.54 


-8.5 


0.54 


34 


3.4 


-1.4 


-10.0 


0.03 


-10.9 


0.05 


-11.1 


0.03 


-12.0 


0.24 


-10.9 


0.08 


-10.0 


0.29 


-8.7 


0.37 




3.4 


-1.4 


-9.9 


0.04 


-11.0 


0.05 


-11.1 


0.04 


-12.1 


0.27 


-11.0 


0.09 


-9.9 


0.36 


-8.9 


0.51 


35* 


-4.4 


-1.7 






















-9.5 


0.57 


-9.2 


0.52 


36* 


-1.8 


-1.8 






















-9.9 


0.46 


-9.5 


0.42 


37* 


3.9 


-2.1 






















-9.7 


0.54 


-9.2 


0.50 


38 


4.3 


-2.1 


-11.0 


0.02 


-10.8 


0.05 


-10.7 


0.03 


-10.7 


0.44 


-10.3 


0.10 


-9.4 


0.38 


-9.0 


0.33 




4.3 


-2.1 


-11.0 


0.02 


-11.0 


0.05 


-10.9 


0.04 


-11.1 


0.44 


-10.6 


0.10 


-10.1 


0.32 


-9.9 


0.28 


39 


-1.5 


-2.5 


-9.6 


0.05 


-9.8 


0.08 


-9.7 


0.06 


-10.1 


0.57 


-10.2 


0.10 


-10.0 


0.28 


-10.0 


0.21 




-1.6 


-2.5 


-9.9 


0.05 


-10.4 


0.08 


-10.4 


0.07 


-10.6 


0.56 


-10.4 


0.12 


-10.1 


0.33 


-10.1 


0.27 


40 


4.2 


-2.5 


-9.9 


0.05 


-10.5 


0.06 


-10.5 


0.04 


-10.4 


0.50 


-10.2 


0.11 


-9.6 


0.34 


-9.1 


0.32 


41 


-1.8 


-2.5 


-9.8 


0.03 


-9.0 


0.11 


-9.1 


0.07 


-9.3 


0.83 


-9.1 


0.17 


-9.3 


0.40 


-9.0 


0.33 




-1.8 


-2.5 


-10.2 


0.04 


-10.3 


0.09 


-10.2 


0.08 


-10.2 


0.65 


-10.0 


0.15 


-9.9 


0.36 


-9.8 


0.32 


42 


-1.5 


-2.6 


-9.6 


0.04 


-9.5 


0.09 


-9.3 


0.07 


-10.0 


0.60 


-8.7 


0.21 


-8.8 


0.48 


-8.4 


0.43 


43 


2.8 


-2.7 


-10.3 


0.03 


-10.4 


0.06 


-10.3 


0.04 


-10.2 


0.56 


-9.9 


0.12 


-8.6 


0.54 


-7.8 


0.59 


44 


1.4 


-2.9 


-9.5 


0.06 


-9.6 


0.14 


-9.4 


0.14 


-9.4 


0.94 


-9.3 


0.22 


-8.8 


0.65 


-8.4 


0.64 


45 


3.3 


-3.0 


-10.1 


0.03 


-9.9 


0.07 


-9.8 


0.05 


-10.1 


0.59 


-9.6 


0.13 


-8.9 


0.47 


-8.5 


0.43 




3.3 


-3.0 


-10.1 


0.03 


-10.1 


0.08 


-10.1 


0.06 


-10.6 


0.55 


-9.8 


0.15 










46 


-1.8 


-3.1 


-9.4 


0.08 


-9.9 


0.08 


-9.9 


0.05 


-9.8 


0.66 


-9.6 


0.13 


-8.6 


0.54 


-8.2 


0.49 


47 


2.8 


-3.2 


-9.9 


0.03 


-9.2 


0.10 


-9.3 


0.07 


-9.8 


0.65 


-8.8 


0.20 


-8.2 


0.64 


-7.7 


0.62 




2.8 


-3.2 


-10.5 


0.03 


-10.2 


0.09 


-10.1 


0.09 


-11.0 


0.46 


-9.6 


0.20 










48 


-0.6 


-3.2 


-10.1 


0.03 


-10.3 


0.09 


-10.1 


0.09 


-10.2 


0.66 


-9.6 


0.20 


-9.2 


0.57 


-9.0 


0.52 




-0.7 


-3.2 


-10.4 


0.03 


-10.6 


0.05 


-10.4 


0.04 


-10.5 


0.47 


-10.1 


0.11 


-9.3 


0.38 


-8.7 


0.38 


49 


0.3 


-3.6 


-10.1 


0.04 


-9.7 


0.08 


-9.6 


0.06 


-9.8 


0.67 


-9.5 


0.14 


-9.6 


0.35 


-9.1 


0.31 




0.2 


-3.5 






















-10.2 


0.37 


-9.1 


0.57 


50* 


1.6 


-3.5 






















-10.5 


0.35 


-10.1 


0.32 


51 


2.7 


-3.6 


-11.6 


0.01 


-11.2 


0.04 


-11.0 


0.03 


-11.2 


0.35 


-10.6 


0.09 


-9.3 


0.39 


-8.1 


0.51 




2.7 


-3.6 


-11.4 


0.01 


-11.2 


0.05 


-11.1 


0.04 


-11.4 


0.38 


-10.5 


0.11 


-9.9 


0.37 


-8.9 


0.57 


52 


2.3 


-3.7 


-9.7 


0.07 


-10.1 


0.12 


-10.1 


0.12 


-10.9 


0.49 


-9.7 


0.21 










53 


0.2 


-3.8 


-9.3 


0.06 


-10.0 


0.07 


-9.9 


0.05 


-10.0 


0.60 


-10.1 


0.11 


-9.9 


0.30 


-9.6 


0.26 




0.1 


-3.8 






















-11.4 


0.20 


-10.9 


0.19 


54 


1.1 


-3.8 


-9.6 


0.06 


-9.9 


0.08 


-9.8 


0.05 


-10.1 


0.57 


-9.9 


0.12 


-9.1 


0.43 


-8.6 


0.40 


55 


1.1 


-4.1 


-9.6 


0.07 


-9.9 


0.08 


-9.9 


0.05 


-11.1 


0.37 


-9.5 


0.15 


-8.8 


0.49 


-8.2 


0.49 



Notes: Candidate star clusters detected with IRAF/DAOPHOT. Col (1) star cluster candidate ID; Col. (2)-(3) projected offset in arcseconds from the photometric nucleus 
(RA 20''37'«14.123^ Dec 66°06'20.09"); Col. (4)-(15) absolute magnitude (m-M=31.9), M^, and eiTor in magnitude, cr^, in the F330W, F547M, F606W, F814W, FllOW and 
F160W band HST images. The candidates with a star in column 1 could not be fitted with an acceptable star cluster model in (see ® . 
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Table 5. Star cluster ages, extinctions and masses 



ID 


Age [Myr] 


E(B-V) 


log(Mass/M0) 




1 


265 


0.44 


5.9 


0.56 


2 


36 


0.02 


4.6 


1.13 


3 


964 


0.00 


6.4 


5.01 


4 


88 


0.38 


5.6 


0.60 


5 


78 


0.35 


5.6 


3.90 


6 


60 


1.15 


5.5 


14.08 


7 


184 


0.80 


6.0 


2.50 


8 


440 


0.62 


6.4 


1.25 


9 


6 


0.31 


4.0 


2.87 


10 


75 


1.14 


5.8 


3.71 


12 


220 


0.72 


6.1 


1.11 


13 


998 


0.93 


6.4 


1.11 


15 


130 


0.90 


6.0 


0.72 


17 


196 


0.50 


5.6 


4.20 


18 


130 


0.23 


5.2 


2.02 


19 


122 


0.40 


5.6 


1.54 


20 


864 


0.00 


6.1 


174.29 


21 


130 


0.56 


5.6 


1.65 


22 


140 


0.41 


6.1 


0.37 


23 


148 


0.83 


5.6 


0.11 


24 


265 


0.00 


6.8 


340.26 


25 


234 


0.65 


5.9 


1.52 


26 


420 


0.43 


5.8 


1.24 


27 


4 


1.26 


4.6 


0.84 


28 


700 


0.49 


6.2 


2.29 


29 


135 


0.30 


5.8 


3.29 


30 


134 


1.00 


6.0 


21.18 


31 


828 


0.78 


6.8 


125.41 


32 


998 


0.75 


6.5 


6.87 


33 


998 


0.44 


6.2 


1.68 


34 


996 


0.90 


6.8 


8.48 


38 


145 


0.51 


5.9 


0.54 


39 


420 


1.04 


6.4 


4.04 


40 


490 


0.60 


6.2 


0.12 


41 


100 


1.04 


5.5 


2.45 


42 


136 


0.31 


5.2 


1.32 


43 


216 


0.42 


5.8 


1.15 


44 


190 


0.69 


5.6 


0.37 


45 


162 


0.50 


5.6 


0.49 


46 


450 


0.36 


5.8 


0.52 


47 


64 


0.64 


5.3 


0.53 


48 


241 


0.21 


5.7 


0.48 


49 


122 


0.43 


5.4 


5.49 


51 


100 


0.18 


5.7 


1.97 


52 


336 


0.72 


6.0 


0.68 


53 


836 


0.62 


6.2 


5.88 


54 


248 


1.04 


6.0 


2.38 


55 


294 


0.58 


5.8 


2.85 



Notes: Results of the;^'^ fitting of age (column 2), color excess (column 3) and mass (column 4) of the observed star clusters. 
Only those candidates which were successfully fitted are listed. The star cluster ID in column 1 is equal to that in Table |4] 
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